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Abstract 

The most general model-independent analysis of the rare B decay, B — > X s l + l~ , 
is presented. There are ten independent local four- Fermi interactions which may 
contribute to this process. The branching ratio, the forward-backward asymmetry, 
and the double differential rate are written as functions of the Wilson coefficients of 
the ten operators. We also study the correlation between the branching ratio and 
the forward-backward asymmetry by changing each coefficient. This procedure tells 
us which types of operator contribute to the process, and it will be very useful to 
pin down new physics systematically, once we have the experimental data with high 
statistics and the deviation from the Standard Model is found. 
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1 Introduction 



Rare l?-meson decays are very useful for constraining new physics beyond the Standard Model 
(SM). In particular, the processes B — > X s ^f and B — > X s l + l~ are experimentally clean, and are 
possibly the most sensitive to the various extensions to the SM because these decays occur only 
through loops in the SM. Nonstandard model effects can manifest themselves in these rare decays 
through the Wilson coefficients, which can have values distinctly different from their Standard 
Model counterparts. (See for example, ]l], [2], ^ [|, ||, |(| [7|, j|].) Compared to B — > X s j, the flavor 
changing leptonic decay B — > X s l + l~ is more sensitive to the actual form of the new interactions 
since we can measure experimentally various kinematical distributions as well as a total rate. 
While new physics will change only the systematically uncertain normalization for B — > X s 7, 
the interplay of various operators will also change the spectra of the decay B — > X s l + l~ . 

We note that the previous studies towards the model independent analysis have been limited 
to the subset of the ten local four- Fermi interactions within the specific extended models, such 
as the two-Higgs-doublet model, the minimal supersymmetric model, the left-right symmetric 
model, and etc. (See for example, ^, [Tcfl.) In the SM and in many of its extensions, the 
decay B — » X s l + l~ is completely determined phenomenologically by the numerical values of 
Wilson coefficients of only three operators evaluated at the scale jx ~ m^. However, the most 
interesting case would be that this three parameter fit is found unsuccessful to explain the real 
experimental distributions, and that the new interactions necessarily involve an extension of the 
full operator basis to include new operators beyond the usual set ||. And, therefore, the new 
physics senario can be much richer than any of those models. 

In the present paper, instead of doing a model by model analysis, we analyse B — > X s l + l~ in 
the most general model-independent way. There are ten types of local four-Fermi interactions 
which may contribute to the process. This is contrasted to the SM in which there are only two 
local interactions coming from the Box, 7 and Z penguin diagrams and one non-local interaction 
coming from the cascade decay, B — > Xgj* — > X s l + l~ [|ll], 12, 13, 14, [HJ. Our final purpose is 
to develop the most general way to distinguish the various four-Fermi interactions, which can 
be fully investigated at B factories of KEK-B, SLAC-B, B-TeV and LHC-B. As is well known, 
all results in direct muon decay (energy spectra, polarizations and angular distributions) and in 
inverse muon decay (the reaction cross section) at energies well below M\\r may be parametrized 



[16] in terms of 10 (complex) couplings, g^n (a, (3 = R,L,j = S,V,T), and the Fermi coupling 
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constant Gp using the matrix element, 

—pr alp < eairVe >< p M |r 7 |/z^ > . (1) 

The "V — A interaction" of the SM corresponds to the single coupling g^ L , which has been 
experimentally probed thoroughly. 

As quantities to be measured in our first step, we study the decay rate B and the forward- 
backward (FB) asymmetry A Jl7|]. These are written as functions of the coefficients of the 
interactions, i.e., £>({Cj}) and A({C{}), where {Cj} is a set of Wilson coefficients. We study 
the branching ratio and the FB asymmetry by changing each coefficient. This shows how 
each operator contributes to the measurements in a different way from the other operators. 
Furthermore, to distinguish types of operators, we draw the correlation between the decay rate 
and the asymmetry. By changing each coefficient, we can draw a flow in the two dimensional 
plane (B(C m ),A(C m )) |l8| . This flow depends on the type of the interaction which acts. Thus 
we may see which type of operator contributes to the process once these quantities are measured 
precisely. 

The paper is organized as follows. In Section 2 we show the most general four- Fermi in- 
teractions and compute the decay distribution. We also derive the branching ratio and the 
asymmetry. Their dependence on the four-Fermi interactions is also studied. In Section 3 the 
correlation between the branching ratio and the asymmetry is discussed and Section 4 is devoted 
to our summary. 

2 The four-Fermi interactions and the decay distributions 

We start by defining the various kinematic variables. In the present paper, the inclusive semilep- 
tonic B decay is modeled by the partonic calculation, i.e., b(pb) — ► s(p s ) + l + {p + ) +Z~(p_). This 
is regarded as the leading order calculation in the l/m& expansion |l9|, pCp . Then the decay 
distribution is described by the following two kinematic variables s and u, 

s = (p b - p s ) 2 = (p + + p-) 2 = ml + m 2 s + m 2 + + m 2 _ - t+ - t_, 

t + = {p s +p + ) 2 = (Pb-P-) 2 , 
t- = (p s +p-) 2 = {pb-p+f, 

u = t + -t— (2) 
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In the center of mass frame of the dileptons, u is written in terms of 6, i.e., the angle between 
the momentum of the B meson and that of l + , 

u = —u(s) ■ cosO = —u(s)z, 
with z = cos 9, 



/ 4l7l? 

u(s) = y(s- (m b + m s ) 2 )(s- (m b -m s ) 2 )(l j-). (3) 

The phase space is defined in terms of s and z, 

4m 2 < s < (m b — m s) 2 , 

-l<z<l. (4) 

The matrix element of the decay is written as the sum of the SM contribution and the 
contribution from the local four- Fermi interactions, 

M=M SM +M NEW , (5) 

where M SM is the SM part and is given by 



M SM = ^v t * s v tb [ (c e 9 ff -c w )s L1 ,b L i L1 n L 

+ (C 9 e// + do) s Llfl b L Tr^Ir 

- 2C e 7 ff sia^(m s L + m b R)b Ij^l ], (6) 

We use the following values for the Wilson coefficients (C^ DR = 4.153, C w = -4.546, C 7 = 
—0.311). They correspond to the next-to leading calculation of QCD corrections |l4|, [l5|]. The 
renormalization point /i and the top quark mass are set to be [i = m b = 4.8 (GeV) and 
rat = 175 (GeV), unless otherwise specified. (We follow Ref. pll] for the choice of the parameters 



in the SM as well as the incorporation of the long-distance effects of charmonium states 22, E 



There are ten independent local four-Fermi interactions which may contribute to the process. 
A1 NEW is a function of the coefficients of local four-Fermi interactions and is defined as 

Gpct - 

a^new = -f=-Kv a [ c LL sli^l i L1 n L 



+ Clr SL^txbL Ir-j^Ir 

+ Crl sr^^r Il^Il 

+ Crr sr^^r Ir^Ir 

+ Clrlr s]J)r IlIr 
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+ Crllr SRbL IlIr 

+ Clrrl sl&r IrIl 

+ Crlrl SRbL IrIl 

+ C T sa^b la^l 

+ iC TE sa^b l<j aP l e^ a % 



(7) 



where the Cx's are the coefficients of four- Fermi interactions. Among them, there are four vector 
type interactions denoted by Cll, Clr, Crl, and Crr. Two of them (Cll, Clr) are already 
present in the SM as the combinations of (Cg — Cio, Cg + Cio). Therefore they correspond to 
the deviation of these coefficients from the SM values. The other interactions, denoted by Crl 
and Crr, are obtained by interchanging the chirality projections L and R. There are four scalar 
type interactions, Clrlr, Crllr, Crllr, and Crlrl- The remaining two, denoted by C T and 
Cte, correspond to tensor type. 

It is now straightforward to compute the double differential decay rate. The results for the 
massless lepton case, B — ► X s e + e~ and B — > X s fi + fi~ , are explicitly written in the text. The 
massive case is separately given in the Appendix. The double differential rate is a function of the 
dilepton invariant mass s and z = cos 9. In the massless limit, the interference between tensor 
(scalar) type interactions and vector type interactions vanishes. Then the expression becomes 
considerably simplified, 



dB 



1 



dsdz 2mb' 



+ A 2 (s,z 

+ A 3 (s,z 

+ At(s,z 

+ A 5 (s,z 



{4|C 7 | 2 } 
{ (C 9 e// - Cm) 
9" ~ Cio) 



+ 



{ (c 9 e// 



(c e 9 ff + c 10 )\ 2 } 

(C e 9 ff + Cm) 2 } 



{2Re(-2C 7 (C e 9 ff * - C* w )) + 2Re(-2C 7 (C e 9 ff * + C* 10 ))} 
{2Re(-2C 7 (C e 9 ff * - C* w )) - 2Re(-2C 7 (C e 9 ff * + C* w ))} 



+ 


A 2 {s, 


z) 


{2Re((C e 9 ff - C l0 )C* LL ) + 2Re((C e 9 ff + C 10 )Cl R )} 


+ 


A 3 (s, 


z) 


{2Re((C e 9 ff - C 10 )C* LL ) - 2Re((C e 9 ff + C w )Cl R )} 


+ 


A A {s, 


z) 


{2Re(-2C 7 (Cl L + C* LIl ))} 


+ 


A 5 {s, 


z) 


{2Re(-2C 7 (Cl L - C* LB ))} 


+ 


A 6 (s, 


z) 


{2 J Re(-2C 7 (C^ L + C3 aR ))} 


+ 


A e {s, 


z) 


{-2Re((C e 9 ff - C 10 )C* RL ) - 2Re((C e 9 ff + C 10 )C^ R )} 
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+ A 7 (s, z) {2Re(-2C 7 (C* RL - C RR ))} 



+ A 2 (s,z) {|C LL | 2 + |C LR | 2 + |C RL | 2 + |C RR | 2 } 

+ A 3 (s, z) {|C LL | 2 - |C LR | 2 - |C RL | 2 + |C RR | 2 } 

+ A 3 (s, z) {-4Re[C LRLR (C^ - 2C^ E ) + C RLRL (C^ + 2C^ E )]} 

+ A 6 (s,z) {-2Re(C LL C RL + C LR C RR ) 

+Re(C LRLR C RLLR + C LRRL C RLRL )} 

+ A 8 (s,z) {|C LRLR | 2 + |C RLLR | 2 + |C LRRL | 2 + |C RLRL | 2 } 

+ A 9 (s,z) {16|C T | 2 +64|C TE | 2 }], (8) 

where A n are the functions of kinematic variables s and z, 

l(m 2 -m 2 ) 2 )}, (9) 
(10) 

(11) 
(12) 

(13) 
(14) 
(15) 
(16) 
(17) 
(18) 

Among them, Ai, A 2 , A^, Aq, Ag, Ag are even functions with respect to z while ^3, A§, and A 7 
are odd functions with respect to z. The former contribute to the single differential rate dB/ds 
and the latter contributes to the forward-backward (FB) asymmetry dA/ds. 
Let us first investigate the single differential rate, 

dB f l d 2 B , 

dz. (19) 



M{s,z) 


= -u(s){-16m 2 m 2 s s + (m 2 + m 2 s )(-2s 2 + 2u(s) 2 z 2 


A 2 (s,z) 


= u (s){-u(s) 2 z 2 -s 2 + (m 2 -m 2 ) 2 }, 


A 3 {s,z) 


= 2u(s) 2 zs, 


A 4 (s,z) 


= u(s){2(m 2 +m 2 ) S -2(m 2 -m 2 ) 2 }, 


A 5 (s,z) 


= —2{m 2 + m 2 )u(s) 2 z, 


A 6 (s,z) 


= 4mbm s su(s), 


A 7 {s,z) 


= 4mbm s u(s) 2 z, 


A 8 (s,z) 


= u(s)(m 2 + m 2 — s)s, 


Ag(s,z) 


= u{s){-2u{s) 2 z 2 - 2s{m 2 b + m 2 ) + 2{m\ - m 2 ) 2 }, 


u(s) 


= ^{s - (m b - m s ) 2 }{s - (m b + m s ) 2 }. 



ds J -1 dsdz 
It is given by 



+ M 2 (a) {2 C 9 e// ' + 2\C 



2 . _. . 2 



'10 1 



6 



+2Re(C e 9 ff (C* LL + C* LR )) + 2Re(C w (-C* LL + C* LR )) 
+ \Cll\ 2 + \C LR \ 2 + \C RL \ 2 + \C RR \ 2 } 

- M 4 (a) {ARe[C* 7 (2C e 9 ff + C LL + C LR )}} 

- M 6 (a) {ARe[C* 7 {C RL + C RR )) 

+2Re(C e 9 ff (C RL + C* RR )) + 2Re(C 10 (-C RL + C* RR )) 
+2Re(C LL C RL + C LR C RR ) - Re{C LRLR C* RLLR + C LRRL C* RLRL )} 
+ M 8 (s) {\C LRLR \ 2 + \C RLLR \ 2 + \C LRRL \ 2 + \C RLRL \ 2 } 
+ M 9 (s){16\C T \ 2 + 64\C TE \ 2 }}, (20) 

where Bo is a normalization factor, 

3q 2 \V t *V tb \ 2 1 
— B s i - - , (21) 

16vr 2 |y cb | 2 f{m c )K(m c ) 



and the phase space factor, /(m c = ^s), and the 0(a s ) QCD correction factor [24], «(m c ), of 
b — > cZf are given by 

/(m c ) = 1 - 8m c 2 + 8m c 6 - m c 8 - 24m c 4 In m c , (22) 

/ - n n 2a s (m fe ) 2 31 2 3 

K(m c ) = 1 ^ — [(tt -— )(l-m c ) + -]. (23) 

For the numerical calculations, we set = 1 and use B s i = 10.4%. M n are functions which 

are obtained by integrating A n (s,z) with respect to z as follows, 

M n (s) = 1^ A n (s,z)dz = 2A n (s,^/l/3), (24) 

where n is 1,2,4,6,8 and 9 and they are explicitly given by 

Afi (a) = 2u{s)-[-lQm 2 b m 2 s + (ml + m 2 )(-2s 2 + -u{sf + 2(m? -m 2 ) 2 )], (25) 
s 3 

M 2 (a) = -2u(s)[^u(s) 2 + s 2 - (m 2 - m 2 ) 2 }, (26) 

M 4 (a) = 2u(s)[2(ml + m 2 )s - 2(ml - m 2 ) 2 ], (27) 

M 6 (a) = 2u(s)[4m b m s s}, (28) 

M B (a) = 2u(s)[s(m 2 1 +m 2 -s)], (29) 

M 9 (a) = 2ii(s)[-^n(s) 2 -2(mg + m 2 )s + 2(mf-m 2 ) 2 ]. (30) 

In order to compare the kinematic functions M n , we plot them in Figure 1. This plot shows 
that Mi is the largest due to its 1/a dependence and that Mq is suppressed by a factor of m s 
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(the strange quark mass). Knowing the kinematical dependence, we plot the differential decay 
rate dB/ds by changing each coefficient Cx in the new physics amplitude -M NEW . We summarize 
how the differential rate dB/ds is changed accordingly. 

• Vector type interactions 

We assume that Cx are real parameters. This means we do not introduce any new physics 
phase in addition to that of the SM given by Arg(V t * Vtb). In Figs. 2,3,4 and 5, we change 
Cll, Clr, Crl, Crr respectively. Among the vector interactions, dB/ds changes the most 
with respect to Cll- This can be explained as follows. By neglecting the term proportional 
to Mq, we can collect the terms coming from Cll, Clr, Crl, and Crr. They are given by 

M2(s)C 2 LL + {M 2 (s)2(C e 9 ff - Cio) + M 4 (s)(-AC 7 )}C LL , (31) 
M2(s)C 2 LR + {M 2 ( s )2(C 9 e// + Cio) + M 4 (s)(-4C 7 )}C LR , (32) 

M 2 (s)C RL , (33) 
M 2 (s)C 2 RR . (34) 

Around s = 5 GeV 2 , Re(C% ,f -C 10 ) ~ 9.5 and Re(C% ,f +C 10 ) ~ 0.4. Then the interference 
term between Cll and Re(Cg^ — Cio) is large and dB/ds is sensitive to the change of 
Cll- The contribution is constructive for Cll = | Cio I an d destructive for Cll = — \Cio\, 
as shown in Figure 1. The effects of the other vector interactions (Clr, Crr, Crl) are 
smaller than that of Cll- To see this dependence clearly, we have also studied the partially 
integrated branching ratio (B) and have shown how it depends on the coefficients. (See 
Figures 8 and 9.) The range for the integration is 1 < s < 8 (GeV 2 ). This range is below 
the J /ip resonance and dB/ds is nearly flat within the range. From these figures, the 
branchingratio depends strongly on Cll while the dependence on the other coefficients is 
rather weak. The contribution of Cll to B is positive (negative) for Cll > (< 0). The 
contribution due to Crl and Crr is positive. 

• Scalar and tensor type interactions 

The scalar and tensor interactions increase dB/ds, as shown in Figures 6 and 7. This 
can be interpreted as follows. As far as we see the sensitivity to only one of the scalar 
operators, i.e., Clrlr,Crllr, Clrrl and Crlrl, the interference term does not act at 
all. And it is enough to study the response of dB/ds to the change of the combinations 
of the coefficients, \C L rlr\ 2 + \Crllr? + \C L rrl\ 2 + \Crlrl\ 2 and 16|C T | 2 + U\C TE \ 2 . 
Because their contribution is always positive, it increases dB/ds. Furthermore, even if two 
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or more scalar interactions act together, the effect of the interference term denoted by 
Mq is suppressed by a factor of m s and can not exceed the positive contribution. Thus 
the major dependence on the coefficients comes through the combination, \Clrlr\ 2 + 
\Crllr\ 2 + \Clrrl\ 2 + \Crlrl\ 2 , even if two or more operators contribute to the process. 
Therefore, dB/ds increases even if more than one scalar and/or tensor type of interaction 
are present. In Figures 8 and 9, we can see that the branching ratio B always increases as 
one of the scalar and tensor interactions act. 
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5 10 15 20 25 

s [GeV 2 ] 

Figure 2 The differential branching ratio for the SM (thick solid line), C'll is — |Cio| (thin line), 
— 0.7 |Cio| (thin dashed line), 0.7 |Cio | (thin dashed line) and |Cio| (thin line). The coefficients of the 
other interactions in AInew are set to zero. 




5 10 15 20 25 

s [GeV 2 ] 

Figure 3 The differential branching ratio for the SM (thick solid line), Clr is — |Cio| (thin line), 
— 0.7 |Cio| (thin dashed line), 0.7 |Ci 1 (thin dashed line) and | C i0 1 (thin line). The coefficients of the 
other interactions in Mnew are set to zero. 
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dB 
ds 



x 10 




10 15 
s [GeV 2 ] 



25 



Figure 4 The differential branching ratio for the SM (thick solid line), Crl is — |Cio| (thin line), 
— 0.7 |Cio| (thin dashed line), 0.7|Cio| (thin dashed line) and |Cio| (thin line). The coefficients of the 
other interactions in AInew are set to zero. 




5 10 15 20 25 

s [GeV 2 ] 



Figure 5 The differential branching ratio for the SM (thick solid line), Crr is — |Cio| (thin line), 
— 0.7 |Ci | (thin dashed line), 0.7 |Ci 1 (thin dashed line) and | C i0 1 (thin line). The coefficients of the 
other interactions in .Mnew are set to zero. 
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5 10 15 20 25 

s [GeV 2 ] 

Figure 7 The differential branching ratio for the SM (thick solid line), 16 |C T | 2 + 64 \C T e\ 2 = 0.5 \C W 
(thin dashed line) and |Cio| 2 (thin solid line). 
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Figure 8 The partially integrated branching ratio B = J 1 ^dsasa lunction of Ct (thick line), Cte 
(thick dashed line) and Cll (dotted dashed line). 
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Figure 9 The partially integrated branching ratio B = j 1 ^ds as a function of Cll (thick solid line), 
Clr (thick dashed line), Crl (thin dashed line), Crr{ thin solid line ) and Clrlr (dotted line). 



13 



Now let us turn to the forward-backward (FB) asymmetry. The normalized asymmetry 
dA/ds is given by 



ds r 1 , d?B , rO , d?B ' 

aS ds Jo dz dslTz + J-l 

where is the FB asymmetry and is expressed as 



(35) 



dA 1 



2mf, 8 



B [ A 3 (s, 1) { (C?/ f - Cio) " - (C 9 e " + Cio) "} 

+ A 5 ( S , 1) {2i?e(-2C 7 (C 9 e// * - C? )) - 2i?e(-2C 7 (C 9 e// * + C? ))} 

+ A 3 (s, 1) {2i?e((C 9 e// - C 10 )C£ L ) - 2Re((C e 9 ff + C W )C* LR )} 
+ Ms, 1) {2Re(-2C 7 (C* LL - C* LR ))} 
+ A 7 (s,l) {2Re(-2C 7 (C* RL -C* RR ))} 

+ M*, i) {I^lI 2 - |C LiJ | 2 - ICmI 2 + |C Rfi | 2 } 

+ A 3 (s, 1) {-ARe[C LRLR {C* T - 2C* TE ) + C RLRL {C* T + 2C^)]}]. (36) 
J 4 ra (s, 1) are the functions which are obtained by integrating A n (s, z) with respect to z, 

A n (s,l) = j A n (s,z)dz - J ^A n (s,z)dz, (37) 

A 3 (s,l) = 2u(s) 2 s, (38) 
A 5 (s,l) = -2{m 2 b + m 2 s )u{s)\ (39) 
A 7 (s,l) = 4m b m s u(s) 2 , (40) 

where n is 3, 5 and 7. These kinematic functions are plotted in Figure 10. Among them, A 7 is 
smaller than the others because it is suppressed by a factor of m s . In the same way as we have 
done for the differential rate, we plot the normalized asymmetry by changing a coefficient of 
vector type interactions within the range —Cio < Cx < C\q. The results are shown in Figures 
11-14. These plots show that the asymmetry is the most sensitive to Cll, which gives a positive 
(negative) contribution for Cll > (< 0). Figures 13 and 14 tell us that C R l (C rr ) decreases 
(increases) the asymmetry. To interpret the results, we collect the terms coming from a vector 
type interaction in the FB asymmetry, 

Ms){CIl + 2(C 9 e// - C W )C LL } - 4A 5 (s)C 7 C L l, (41) 



yeff 



14 



Ms){-C 2 LR - 2(C 9 e// + C 10 )C LR } + 4A 5 (s)C 7 C L r, 

Ms){-c 2 RL }, 
Ms){c 2 RR }, 



(42) 



(43) 



(44) 



where we neglect the term proportional to Aj. Because (Cg — Cio) is much larger than the 



other SM coefficients, the interference term with Cll is important. Therefore the asymmetry 
is sensitive to Cll- The response to Clr is more involved and the interference term between 
CV and Clr seems to be the most important. It gives a positive (negative) contribution to 
the asymmetry for Clr > (< 0). We also note that Crl (Crr) gives a negative (positive) 
contribution to the asymmetry and it is consistent with the results shown in Figures 13 and 14. 

We now study the integrated FB asymmetry (A). We choose the same integration region 
as was chosen for the branching ratio (B), and define the integrated FB asymmetry A and the 
normalized FB asymmetry A, 



where B = J® ^ds. We outline how A and A depend on the interactions. In Figure 15 (16), 
the dependence of the integrated asymmetry A (normalized asymmetry A ) on \Cx\ is shown 
by changing the coefficient \Cx\ within the small range, — C\o < \Cx\ < Cio- Within this range, 
the linear dependence on the coefficients is important. At first glance, we see that only Cll 
can give rise to both large positive and negative contributions. Clr also gives rise to both 
positive and negative contributions. Crr gives a positive contribution and Crl gives a negative 
contribution as already discussed for dA/ds. We also note that the integrated FB asymmetry 
does not change at all as we change a coefficient of one of the scalar or tensor interactions. (See 
Figure 15.) This is because the dependence on these types of interactions comes only from the 
interference term. On the other hand, the magnitude of the averaged FB asymmetry is reduced 
slightly if one of these interactions act. The reason is that the denominator of the averaged FB 
asymmetry, B, always increases as scalar and tensor interactions are switched on. 




(45) 



(46) 
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Figure 10 A 3 (s, 1) (thin dotted line), — A 5 (s, 1) (thin solid line ) and A 7 (s, 1) (thick solid line). 
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Figure 12 The normalized FB asymmetry for the SM (thick solid line), Clr is — |Cio| (thin line) and 
| Cio | (thin dashed line). The coefficients of the other interactions in .Mnew a re set to zero. 
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Figure 16 The contribution to the integrated normalized FB asymmetry A from Cll (thick solid line), 
Clr (thick dashed line), Crl (thin dashed line), Crr (thin solid line) and Clrlr (dotted line). 

3 Correlation between the branching ratio and the forward- 
backward asymmetry 

In the previous Sections we have studied how the branching ratios and the asymmetries depend 
on each interaction. In order to distinguish which interaction contributes to the process, we plot 
the correlation between the branching ratio and the asymmetry. By changing a coefficient, we 
can draw a flow in the two dimensional planes (£>, A) and (£>, A). The flows depend on the type 
of interaction which acts. Then we may pin down the type of interaction which contributes to 
the processes once these quantities are measured. We allow the coefficients to vary in a wider 
range, i.e., —20 < Cx < 20. The correlation plots are shown for (-6,-4) in Figure 17, and for 
(B, A) in Figure 18. B, A, and A are also plotted as functions of the various coefficients in the 
range. (See Figures 19-21.) There are uncertainties in the SM predictions which come from the 
renormalization point (/x) dependence of QCD correction. The predictions of the SM are shown 
for three values of thr renormalization point. The flows in the (B,A) plane near the point of 
the SM prediction is summarized as follows. 

• The figures show that the flow for Cll is distinctive from the others. In particular, if Cll 
is negative, the branching ratio decreases substantially while the other interactions tend 
to increase the branching ratio. 
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• As \Clr\ and \Crl\ increase, the asymmetry decreases except the range < Clr < 3. 
The branching ratio tends do increases as \Clr\ and \Crl\ are larger. 

• As \Crr\ increases, the asymmetry and the branching ratio increase. 

• The tensor and scalar interactions do not affect the asymmetry and increase the branching 
ratio. Therefore the flow becomes flat. 

• The branching ratio and the asymmetry does not strongly depend on the sign of Crl and 
Crr. Therefore, the flow for the positive coefficent and that of the negative coefficent are 
nearly degenerate. The flows for the positive Cll and Clr are distinct from those for 
the negative coefficients. They are not degenerate because the branching ratio and/or the 
asymmetry depend on their sign. 

The flow in the plane (B,A) is similar to that in (B,A) for small changes of the coefficients. 
As coefficients are increased, nonlinear dependence begins to be important. This causes the 
difference between A and A. As an example, increasing the coefficient of a scalar or tensor 
operator, in the (B, A) plane, the asymptote of the flow is A = 0, while in the (B,A) plane, it 
is flat. 
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Figure 17 The flows in (B,A) plane. In each flow, Cll (thick solid line ), Clr (thick dashed line), Crl 
(thin dashed line) , Crr (thin solid line) , and the tensor and scalar interactions (dotted line) are changed 
respectively. The predictions of the SM are denoted by A, O, and vy corresponding to three values of the 
renormalization point (yu). 
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Figure 18 The flows in (B,A) plane. In each flow, Cll (thick solid line), Clr (thick dashed line), Crl 
(thin dashed line), Crr (thin solid line), and the tensor and scalar interactions (dotted line) are changed 
respectively. The predictions of the SM are denoted as A, O, and y corresponding to three values of the 
renormalization point (fj,). 
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Figure 21 The contribution to the integrated normalized asymmetry A from Cll (thick solid line), Clr 
(thick dashed line), Crl (thin dashed line), Crr (thin solid line) and the tensor and scalar interactions 
(dotted line). 
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4 Summary 



We have given the most general model-independent analysis of the rare B decay process B — ► 
X s l + l~ . This process is experimentally very clean, and possibly most sensitive to the various 
extensions of the Standard Model (SM), compared to other rare decay processes. As is well 
known, the main reason for studying rare B decays is to measure the effective FCNC vertices in 
order to test the SM precisely and to search for new physics beyond the SM. We have investigated 
the decay B — * X s l + l~ in the full operator basis, instead of doing a model by model analysis. The 
sensitivity to the coefficients of ten independent local four-Fermi interactions is systematically 
studied for the branching ratios and the forward-backward (FB) asymmetries. The correlation 
between the branching ratio and the FB asymmetry has been studied, and the flow for Cll, 
the coefficient of the operator (sL^^bi /l7^l)> is found to be distinctive from the other vector 
type interactions. The reason for the strong sensitivity to Cll comes from the large interference 
between (Cg — Gxo) and Cll- The tensor and scalar operators increase the breaching ratio, while 
the asymmetry is not changed at all. The other flows have also been discussed in detail. 

As further comments, the uncertainties from higher-order QCD corrections can be reduced 
and such a reduction is very useful to distinguish the effect of new physics from the uncertainties 
in the predictions of the SM. Unless the theoretical predictions for the rates and the asymmetries 
in the SM are on firm ground, it would be difficult to make sure that new physics has indeed been 



discovered. A complete short distance NLO calculation for B — > X s l + l is available [14 



and the l/?n| corrections have also been known [17, [H| for some time, but these are found 
to give only very small corrections. The remaining uncertainties are the non-perturbative long 



distance contributions associated with the J/ip and ip' resonances [22|, where some modeling 
uncertainties remain, and the related 1/m^ corrections. 

The other useful measurements, such as the tau polarization asymmetry, CP-violating rate 



asymmetry [23| and triple momentum-spin correlations can be studied and independent infor- 



mation on the coefficients can also be obtained [25|. The analysis based on a specific model may 



be systematically organized [^] in our scheme, because the comparison of the predictions of the 
various models will be easy in our framework. 
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A Appendix 



In this Appendix, we derive the various distributions for the massive lepton case. The formulae 
can be applied to the process B — > X s t + t~ . 

The double differential rate is a function of the dilepton invariant mass s and z = cos6, 



dB 



1 



dsdz 2m\f 



r6n 



[ M 

+ A 2 

+ ^3 

+ A A 

+ A 5 

+ Li 

+ A 2 

+ ^3 

+ A 4 

+ A 

+ A 6 

+ A 6 

+ A 7 

+ L 2 

+ L 2 

+ L 3 

+ L 4 

+ L t 

+ L 5 

+ L 5 

+ L 6 



s, z 

s, z 
s, z 
s, z 
s, z 
s, z 

s, z 
s, z 
s, z 
s, z 
s, z 
s, z 
s, z 
s, z 
s, z 
s, z 
s, z 
s, z 
s, z 
s, z 
s, z 



{4|C 7 | 2 } 

{|(c 9 e// -c 10 )| 2 + |(c 9 e// + c 10 )| 2 } 
{|(c 9 e// -c 10 )| 2 -|(c 9 e// + c 10 )| 2 } 

{2Re(-2C 7 (C e 9 ff * - CJo)) + 2Re(-2C 7 (C e / f * + C* 10 ))} 
{2Re(-2C 7 (C e 9 ff * - a )) - 2Re(-2C 7 (C e 9 ff * + C* 10 ))} 
{2Re((C e 9 ff -C 10 )(C e 9 ff * + C* 10 ))} 

{2Re((C e / f - C 10 )C* LL ) + 2Re((C e 9 ff + C 10 )C* LR )} 
{2Re((C e 9 ff - C 10 )Cl L ) - 2Re((C e 9 ff + C w )Cl K )} 
{2i?e(-2C 7 (C3 LL + C3 LR ))} 
{21te(-2C 7 (Cj- L -Cj- R ))} 
{2i?e(-2C 7 (C3 aL + C^ R ))} 

{-2Re((C e 9 ff - C 10 )C RL ) - 2Re({C e 9 ff + C 10 )C RR )} 



{2Re(-2C 7 (C* T 



RL 



RR 



))} 



{Am b Re(-2C 7 (Cl KLR + C* LKRL )} 
{4m s i?e(-2C 7 (C RLLR + C RLRL )} 
{2Re(-2C 7 {C* T )} 
{2Re(-2C 7 (C* TE )} 

{2Re((C e 9 ff - C 10 )CL R ) + 2Re((C e 9 ff + C 10 )C* LL )} 
{-2Re((C e 9 ff - C 10 )C RL ) - 2Re((C e 9 ff + C 10 )C RR )} 
C 10 )C RR ) + 2Re((C e 9 ff + C 10 )C RL )} 



{2Re((C e /f 

ieff 



LRRL 



{2Re({C e 9 " -Cio)(Cj, RLR 
-2i?e((C 9 e// + Ci )(CJ 
+ L 2 (s,z) (- m6 ){2 J R e ((C 9 e// -C 10 )(Ci RLR + Ci RRL ) 



LRLR ~ C LRRL)} 



+2^e((C 9 e// + Ci )(C 



LRLR 



+ c 



LRRL 



)} 
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+ L 7 (s, z) {2Re{{C e 9 ii - C 10 )(C RLRL - C RLLR ) 
-2Re((C e 9 ff + C 10 )(C RLRL - C RLLR )} 
+ L 2 (s,z) (-m s ){2i?e((C 9 e// -C 10 )(C RLLR + C RLRL ) 
+2Re((C e 9 ff + C 10 )(C RLLR + C RLRL )} 

+ L 8 (s, z) {2Re((C e 9 ff - C W )C* T ) + 2Re((C e 9 ff + C W )C* T )} 

+ L 2 (s, z) I2(m b + m s ){2Re{(C e 9 ff - C 10 )C* T ) - 2Re{(C e 9 ff + C 10 )C* T )} 

+ L 9 (s, z) {2Re({C e 9 ff - C W )C* TB ) + 2Re((C e 9 ff + C 10 )C^ E )} 

+ L 2 (s, z) 2A{-m b + m s ){2Re((C e 9 ff - C 10 )C* TB ) - 2Re((C e 9 ff + C 10 )C^ E )} 

+ A 2 (s,z) {|C LL | 2 + |C LR | 2 + |C RL | 2 + |C RR | 2 } 
+ A 3 (s,z) {|C LL | 2 - |C LR | 2 - |C RL | 2 + |C RR | 2 } 
+ M*, *) {-2Re(C LL C RL + C LR C RR ) 

+Re(C LRLR C RLLR + C LRRL C RLRL )} 

2 2 2 2 

+ M s , z ) {I c lrlrI +I c rllrI + I c lrrlI +I c rlrlI } 

+ A 3 (s,z) {-4Re(C L R LR )(C^-2C^ E )-4Re(CR L R L )(C T + 2C^ E )} 

+ A 9 (s,z) {16|C T | 2 } 

+ A 9 (s,z) {64|C TE | 2 } 

+ Li(a, z) {2Re(C LL Ci R + C RL C RR ) 

-Re(C LRLR CL RRL + C RLLR C RLRL )} 
+ Ms, z) {-2Re(C LL C RL + C LR C RR ) 

+Re(C LRLR C RLLR + C LRRL C RLRL )} 
+ Ms, z) {2Re(C LL C RR + C LR C RL ) 

+-Re(C LRLR C RLRL + C LRRL C RLLR )} 
+ Ms, z) {2Re(C LL - C LR )(C£ RLR - C£ RRL ) 

+2Re(C RL - C RR )(C RLLR - C RLRL )} 
+ Ms, z) (-m b ){2Re(C LL + C LR )(C£ RLR + C£ RRL ) 
+2Re(C RL + C RR )(C 

RLLR + C RLRl)} 
+ Ms, z) {2Re(C LL - C LR )(C RLRL + 

c rllr) 

+2Re(C RL - Crr)(Clrr L - c Lrlr)} 
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+ L 2 (s,z) (-m s ){2Re(C LL + C LR )(C| lLRL + C^ LLR ) 
+2Re(C RL + C rr )(Cl RRL - Cl RLR )} 
+ L 8 (s,z) {2Re(C LL + C LR )(C^) + 2Re(C RL + C RR )(C^)} 
+ L 2 (s, z) 12(m 6 + m s ){2Re(C LL - C LR )(C^) + 2Re(C RL - C RR )(C^)} 
+ L 9 (s, z) {2Re(C LL + C LR )(C^ E ) + 2Re(C RL + C RR )(C^ E )} 
+ L 2 (s, z) 24{-m b + m s ){2Re(C LL - C LR )(C^ E ) 

+2Re(C RL -C RR )(C^ E )} 
+ L 5 (s,z) {-192|C TE | 2 }], (A-l) 



u(s) = y (s - (m b + m s ) 2 )(s - (m b - m s ) 2 )(l -*-), 

A\ = —u(s)[—\Qmlm 2 s(2m 2 + s) 

+K + m 2 s ){mf (-8(m 2 + m 2 s )s + 8(m 2 b - m 2 s ) 2 ) + 2s(-s 2 + u{s) 2 z 2 + (m 2 - m 2 ) 2 )}], 

A 2 = u ( s )(-u(s) 2 z 2 -s 2 + (m 2 -m 2 ) 2 ), 

A 3 = 2u(s) 2 zs, 

A A = u(s)^{2(2m 2 + s)((m 2 + m 2 ) S -(m 2 -m 2 ) 2 ), 

A 5 = -2(m 2 b +m 2 )u(s) 2 z, 

Aq = 4u(s)m b m s (2m 2 + s), 

Ay = 4m b m s u(s) 2 z, 

A 8 = -u(s)(m b +m 2 - s)(2m 2 - s), 

Ag = u(s){4m 2 (m 2 — Qm b m s + m 2 — s) 

-2u{s) 2 z 2 - 2{m 2 b + m 2 )s + 2{m 2 b - m 2 ) 2 }, 

L\ = 4u(s)m 2 (ml + m 2 — s) , 

L 2 = miu(s) 2 z, 

L3 = u(s)-mi(m b + m s )(8s 2 + 8s(m 2 + 6m b m s + m 2 ) — 16(m 2 — m 2 ) 2 ), 

L4 = u(s)-mi(m b — m s )(— 16s 2 — 16s(m 2 — 6m b m s + m 2 ) + 32(m 2 — m 2 ) 2 ), 

L 5 = -I6u{s)m b m s ml 

Lq = u(s)mim b (s - ml + m 2 ) , 

L 7 = u(s)mim s (s + ml - m 2 s ), 



where A n and L n are functions of the kinematic variables s and z, 
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L 8 = 12u(s)mi(-(m b + m s )s + (m b - m s )(ml - m 2 s ), 

Lg = 24u(s)mi((mb — m s )s — (mi> + m s )(m 2 — m 2 s ). (A-2) 

Among them, A±, A 2 , A4, Aq, As, Ag, L\, L3, L4, L5, L$, L7, L% and Lg are even functions with 
respect to z, while A%, A^, A 7 ,L 2 are odd functions. In the limit of mi — 0, L n vanish and A n 
lead to those in the massless limit. 

To investigate the differential decay rate and the FB asymmetry, we define the following 
functions. They are obtained by integrateing the A n (s, z) (n is 1 to 9) and L n (s, z) with respect 
to z as follows, 

Mi{s) = J 1 ^A i {s,z)dz = 2A i (s,Jl/i), 

Nj(s) = J Lj(s,z)dz = 2Lj(s,^Jl/3), (A-3) 
where i is 1,2,4,6,8 and 9, and j is 1,3,4,5,6,7,8 and 9. Also 

A k (s,l) = A k (s,z)dz - J ^A k {s,z)dz, 

£2(s,l) = J L 2 (s,z)dz- J L 2 (s,z)dz, (A-4) 

where k is 3,5 and 7. 

Finally, the differential decay rate is given by 

<iB 1 r£o [ Mx( s ){4|CV| 2 } 



ds 2nib 



+ M 2 (s) { (C e 9 ff - C10) " + (C e / J + C10) "} 

+ M 4 (s) {2Re(-2C 7 (C e 9 ff * - C* 10 )) + 2Re(-2C 7 (C e 9 ff * + C? ))} 

+ MM {2i?e((C 9 e// - C 10 )(C 9 e// * + CJo))} 

+ M 2 ( S ) {2i?e((C 9 e// - C 10 )CL L ) + 2i?e((C 9 e// + C 10 )C£ R )} 

+ M 4 ( S ){2i?e(-2C 7 (Cj jL + Cj jR ))} 

+ M 6 ( S ) {2i?e(-2C 7 (C RL + C RR ))} 

+ M 6 ( S ) {-2 J Re((C 9 e// - C 10 )C RL ) - 2Re((C e 9 ff + C 10 )C RR )} 

+ N 3 {s) {2Re(-2C 7 (C* T )} 

+ N 4 (s) {2Re(-2C 7 (C* TE )} 

+ iVi(s) {2 J Re((C 9 e// - C 10 )C£ R ) + 2i?e((C 9 e// + C W )C* LL )} 

+ N 5 (s) {-2Re((C e 9 ff - C 10 )C RL ) - 2Re((C e 9 ff + C 10 )C RR )} 
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2 



+ N 5 (s) {2Re((Cfl - C 10 )C RR ) + 2Re((Cf f + Cio)C RL )} 
+ N 6 (s) {2Re((C e 9 ff - C 10 )(C£ RLR - C£ RRL ) 

-2Re((C e 9 ff + C 10 )(C£ RLR - Cl RR l)} 
+ N 7 (s) {2Re((C e / f - C 10 )(C RLRL - C RLLR ) 

-2Re((C e 9 ff + C 10 )(C RLRL - C RLLR )} 
+ N 8 (s) {2Re((Cl ff - C 10 )C* T ) + 2Re({C e 9 n + C W )C* T )} 
+ N 9 (s) {2Re((C e 9 ff - C 10 )C^ E ) + 2Re((C e / f + C 10 )C^ E )} 

+ M 2 (s) {|C LL | 2 + |C LR | 2 + |C RL | 2 + |C RR | 2 } 
+ M 6 (s) {-2Re(C LL C RL + C LR C RR ) 

+Re(C LRL RC RLLR + C L rrlC rlrl )} 
+ M 8 (s) {|C LRLR | Z + |C RLLR | Z + | 

c lrrlI + I c rlrlI } 

+ M 9 (s) {16|C T | 2 } 

+ M 9 (s) {64|C TE | 2 } 

+ iVi( S ){2Re(C LL C£ R + C RL C RR ) 

-Re(C LRLR CL RRL + C RLLR C RLRL )} 
+ N 5 (s) {-2Re(C LL C RL + C LR C RR ) 

+Re(C L RLRC RLLR + C LRRL C RLRL )} 
+ iV 5 (s){2Re(C LL C RR + C LR C RL ) 

+ -Re(C LRLR C RLRL + C LRRL C RLLR )} 
+ N 6 (s) {2Re(C LL - C LR )(C£ RLR - C^ RRL ) 

+2Re(C RL - C RR )(C RLLR - C RLRL )} 
+ N 7 (s) {2Re(C LL - C LR )(C RLRL + C RLLR ) 

+2Re(C RL - C rr )(Cl RRL - Cl RLR )} 
+ iV 8 ( S ) {2Re(C LL + C LR )(C^) + 2Re(C RL + C RR )(C^)} 
+ N 9 (s) {2Re(C LL + C LR )(C^ E ) + 2Re(C RL + C RR )(C^ E )} 
+ N 5 (s) {-192|C TE | 2 }]. (A-5) 
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And the FB asymmetry is 
dA 1 



ds 2m b 8 



[ A 3 
+ A 5 



+ 


As 


(M) 
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A 5 


(M) 
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A 7 




+ 


L 2 
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L 2 


5,1) 
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L 2 





+ L 2 

+ L 2 
+ L 2 



+ 


As 




+ 


As 




+ 


L 2 ( 


5,1) 



+ L 2 

+ L 2 
+ L 2 



S , 1) { (C e 9 ff - C 



10, 



(c 9 e// + c 10 )-} 



S , 1) {2i?e(-2C 7 (C 9 e// * - C? )) - 2 J Re(-2C 7 (C 9 e// * + C? ))} 



1* 1* 

RL ~ °RRi 



))} 

LRLR + C LRRL)} 
RLLR + C RLRL)} 

f - Cio)(Clrlr + c lrrl) 
+2Re((C e 9 ff + C 10 )(C3 LRLR + CL RRL )} 
s, 1) (-m s ){2Re((C e 9 ff - C 10 )(C RLLR + C RLRL ) 
+2Re((C e 9 ff + C 10 )(C RLLR + C RLRL )} 

9 - C 10 )C* T ) - 2Re((C e 9 ff + C 10 )C* T )} 



s, 1) 12(m b + m s ){2Re{{C 9 }f 
s, 1) 24(-m b + m s ){2Re((C e 9 ff - C 10 )C* TE ) - 2Re((C e 9 ff + C 10 )C^ E )} 



{|c L lI 2 -|c L rI 2 -|c rl | 2 + |c rr | 2 } 

{-4Re(C LRLR )(C^ - 2C^ E ) - 4Re(C RLRL )(C T + 2C^ E )} 
(-m 6 ){2Re(C LL + C LR )(CL RLR + C^) 

+2Re(C RL + C RR )(C RLLR + C RLRL )} 
(-m s ){2Re(C LL + C LR )(C RLLR + C RLRL ) 

+2Re(C RL + C RR )(C£ RLR + C£ RRL )} 
12(m b + m s ){2Re(C LL - C LK )(C* T ) + 2Re(C RL - C RR )(C^)} 
24(-m b + m s ){2Re(C LL - C LR )(C^ E ) 

+2Re(C RL - C RR )(C^ E )}]. (A-6) 



s,l) 

8,1) 
8,1) 
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